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ABSTRACT 

We report on the properties of 71 known cataclysmic variables (CVs) in photomet- 
ric Ha emission line surveys. Our study is motivated by the fact that the Isaac Newton 
Telescope (INT) Photometric Ha Survey of the northern galactic plane (IPHAS) will 
soon provide r', i' and narrow-band Ha measurements down to r' ~ 20 for all northern 
objects between —5° < b < +5°. IPHAS thus provides a unique resource, both for 
studying the emission line properties of known CVs and for constructing a new CV 
sample selected solely on the basis of Ha excess. Our goal here is to carry out the first 
task and prepare the way for the second. In order to achieve this, we analyze data on 
19 CVs already contained in the IPHAS data base and supplement this with identical 
observations of 52 CVs outside the galactic plane. 

Our key results are as follows: (i) the recovery rate of known CVs as Ha emitters 
in a survey like IPHAS is ~ 70 per cent; (ii) of the ~ 30 per cent of CVs which were 
not recovered ~ 75 per cent were clearly detected but did not exhibit a significant Ha 
excess at the time of our observations; (iii) the recovery rate depends only weakly on 
CV type; (iv) the recovery rate depends only weakly on orbital period; (v) short-period 
dwarf novae tend to have the strongest Ha lines. These results imply that photometric 
emission line searches provide an efficient way of constructing CV samples that are 
not biased against detection of intrinsically faint, short-period systems. 

Key words: surveys - binaries: close - novae, cataclysmic variables 



1 INTRODUCTION 

Cataclysmic variables (CVs) are interacting binary systems 
in which a white dwarf (WD) primary accretes matter from 
a main-sequence secondary via Roche lobe overflow. If the 
magnetic field of the WD is dynamically unimportant, the 
accretion process takes place entirely via a disc surrounding 
the WD. By contrast, if the magnetic field is very strong, the 
accretion stream from the secondary is channelled directly 
onto the magnetic poles of the WD. Finally, in the interme- 



diate range of field strengths, a partial disc may form that 
is truncated on the inside by the WD magnetic field. 

An observational feature that is common to all of these 
different accretion modes are (Balmer) emission lines. In 
weakly-magnetic CVs, the observed line emission may orig- 
inate frranoptically thin or irradiated parts of the accretion 
disc iWilliamall98Cl) . In magnetic CVs, line emission may 
be produced in the accretion stream, any residual accre- 
tion disc, or in the so-called accretion curtains that chan- 
nel material from the inner edge of a truncated disc to the 
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magnetic poles of the WD. The irradiated secondary star 
in CVs can also contribute an appreciable contrib ution to 
the ob served line emission, particularly in polars. IWarnerl 
( 1995) provides a comprehensive review of CVs, including 
their spectroscopic properties. 

Given the ubiquity of line emission amongst CVs, emis- 
sion line surveys offer a powerful way to fin d new CVs (for 
a cur ren t example of such a CV se arch, see lGansicke et ail 
l2002l and lAungweroiwit et al.l2005Tl . What makes this strat- 
egy particularly promising is that, empirically, the intrin- 
sically faintest, low mass transfer rate (M) systems tend 
to have the lar gest Balmer line equivalent widths (EWs; 
|Pattersonlll984h . Population synthesis models suggest that 
such low accretion rate systems should totally dominate 
the galactic CV population and should be found predomi- 
nantly at short orbital periods, that is below the well-known 
CV "period gap" between 2.2 hrs and 2.8 hrs (|Kolblll993t 
iHowell et alJll99^ 1. However, this dominant population of 
faint, short-period CVs has proven quite elusive, and it is 
still not clear whether this is due to selection effects or 
whether it points to a serious flaw in our understanding 
of CV evolution. Since emission line surveys should be very 
good at finding low M CVs, they should be an excellent way 
of detecting the large population of faint, short-period CVs, 
if it exists. 

The INT/WFC Photometric Ha Survey of the north- 
ern galactic plane (IPHAS) is currently surveying the Milky 
Way in r' ,i' and Ha and provides an excellent data base for 
a detailed CV search at low galactic latitudes. The survey 
goes to a depth of r' ~ 20 and covers the latitude range 
—5° < b < +5°. A detaile d introduction to the survey is 
given bv lDrew etaH l|2005h . 

In this paper, we investigate the properties of the known 
population of CVs that have been observed by IPHAS, sup- 
plemented with additional observations of CVs outside the 
galactic plane with the IPHAS set-up. Such a study is use- 
ful for two main reasons. First, it allows us to determine the 
recovery rate of CVs as Ha emitters as a function of CV sub- 
class, orbital period and apparent magnitude. These recov- 
ery rates effectively specify the completeness of the survey 
with respect to CVs, which is crucial for the interpretation of 
the new CV sample that is currently being constructed based 
on IPHAS data. Second, the existence of a large, uniform 
sample of broad- and narrow-band magnitudes for known 
CVs allows us to check for and update correlations between 
key observables (for example Ha excess) and intrinsic CV 
parameters (for example orbital period and absolute magni- 
tude). 

The structure of this paper is as follows. Section |5] dis- 
cusses the observations obtained. Section [£] discusses the 
two CV samples analysed in the paper and how they have 
been observed. The techniques for selecting CVs from pho- 
tometric data are presented in section [4] In section [£] we 
present and discuss our estimates for CV recovery /detection 
rates with IPHAS-like emission line surveys. In section |5] 
we consider correlations between key observables and sys- 
tem parameters. Our findings are discussed in section |7| 
focusing particularly on the ability of IPHAS to uncover the 
long-sought population of faint, short-period CVs. We sum- 
marize our conclusions in section |H| 



2 OBSERVATIONS 

Our study relies on two observational data sets. The first 
is composed of the IPHAS galactic plane observations that 
were completed up to and including July 2004; this amounts 
to about 1/3 of the final survey area. The second consists 
of a series of imaging observations obtained in June 2004 of 
northern CVs with galactic latitudes |6| > 5°. Both datasets 
were obtained with the same instrumental set-up and ob- 
serving strategy, and both were reduced with the standard 
IPHAS pipeline. 

Full details of the IPHAS ob serving strategy an d reduc- 
tion procedures may be found in lDrew et al.l J2005f) . Briefly, 
all observations were obtained using the Wide Field Cam- 
era (WFC) on the Isaac Newton Telescope, which gives a 
spatial pixel size of 0.333"x 0.333" over a field of view of 
approximately 0.3 square degrees. Photometry was carried 
out using a set of three filters, comprising a narrowband Ha 
filter and an additional set of two broadband Sloan r' and i' 
filters. Exposure times were 120s for the Ha images, 10s for 
the i! band images, and 10s [30s] seconds for all r' band im- 
ages obtained before [after] June 2004. Finding charts from 
iDownes et all fcOOll) were used to identify the CVs in our 
images. 

The final IPHAS data base will contain only observa- 
tions that meet a strict set of quality constraints. However, 
in order to avoid restricting our CV sample too much in the 
present analysis, we did not limit ourselves to fields that have 
passed the IPHAS quality controls. Instead, for each CV, we 
simply took the data from the best available set of obser- 
vations. 1 . Thus some of the data analysed here will not be 
included in the final published survey catalogues. Instances 
in which the data quality of a field was below the nominal 
limits will be indicated explicitly. However, our results from 
these fields are consistent with the results from fields with 
data quality that passes the IPHAS quality controls. 

We finally note that, in all of our analysis below, we 
only work with sources brighter than a pre-defined magni- 
tude limit of r' — 19.5. In principle, IPHAS goes somewhat 
deeper than this, to about r' — 20. However, the increased 
photometric scatter exhibited by the very faintest sources 
makes it counter-productive to include these sources in our 
selection algorithm (see Section 0J. 



3 CONSTRUCTION OF CV SAMPLES 

As noted above, the data we use in our study come from two 
distinct sources. The first is the standard IPHAS database, 
the second is a set of observations of CVs at higher galac- 
tic latitudes observed with the IPHAS set-up. The majority 
of the fields in the second data set exhibit a significantly 
lower object density than the IPHAS fields. This may be 
expected to affect CV recovery rates, and we therefore re- 
tain the distinction between In-Plane and Off-Plane samples 
in our study. 



1 Note that a given CV can be present in multiple IPHAS point- 
ings because (i) most locations are observed at least twice with 
the IPHAS tiling pattern, and (ii) poor quality observations may 
already have been repeated. 
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3.1 The In-Plane CV Sample 

Adopting a limiting magnitude of r' = 19.5 the IPHAS 
imaging database as of July 2004 contains 19 known cat- 
aclysmic variables 2 l isted in the Ritter & Kolb catalogue 
jRitter fc Kolbll2003ft 3 . We extracted standard IPHAS cat- 
alogues, containing object positions and magnitudes, for all 
of the IPHAS fields that include these CVs. As noted above, 
when multiple observations were available for a given CV, 
we always worked with data from the field observed under 
the best conditions. Information and data on each CV in the 
sample are listed in Table [5] This includes a flag to mark 
data that were extracted from fields that do not meet the 
nominal IPHAS quality thresholds. 

3.2 The Off-Plane CV Sample 

The targets for off-ga lactic-plane (\b \ > 5 °) observations 
were again chosen from lRitter fc Kolbl J2003), with the main 
additional constraint being that the target should be observ- 
able in June 2004. A sample of 52 Off-Plane CVs turned out 
to satisfy our magnitude limit of r < 19.5. The IPHAS 
pipeline-reduced data were extracted in exactly the same 
way as for the In-Plane CVs. 

One target, the dwarf nova (DN) IY UMa, deserves ad- 
ditional comment. It was found to have a faint close com- 
panion which is not identified in the finding charts and has 
not been previously described in literature. In the r' band, 
the faint companion has a magnitude of 19.27 mag and is 
separated from IY UMa by 3.5 arcsec. The faint compan- 
ion has resulted in IY UMa being detected as an unresolved 
blend in the i' band, which makes the il band magnitude 
unreliable in this case. 



4 SELECTING Ha EXCESS OBJECTS 

We have developed a systematic selection algorithm to iden- 
tify Ha excess outliers (and hence likely emission line ob- 
jects) from r' — Ha vs r' — i! colour-colour plots of in- 
dividual fields. The selection process does not require user 
intervention and allows us to separate the CVs into those 
that show clear Ha emission and those that do not. We can 
then use this information to test for statistical completeness 
and selection biases in samples of CVs constructed solely on 
the basis of photometric Ha excess. In the following subsec- 
tions, we describe the sequence of steps used by our selection 
algorithm to identify likely emission line sources. 

4.1 Initial Selection Cuts 

The first step in our selection procedure is to apply the fol- 
lowing criteria to all objects in the catalog for each field 
in the IPHAS database. First, objects must be detected in 
all three photometric bands and must not fall on bad pix- 
els of the CCDs. Second, the objects must have positions in 

2 Note that we omit the double-degenerate AM CVn systems 
from consideration; these He-rich systems belong to a different 
evolutionary channel and are not expected to show Ha emission. 

3 The Ritter & Kolb catalogue is available online at 
|http://physics.open.ac.uk/RKcat| 



the three photometric bands that match to within 1 arcsec. 
Third, objects must be classified as "stellar" in the i' band, 
and "stellar" or "probably stellar" in the r' and Ha bands. 
We allow "probably stellar" classification in r' and Ha to 
avoid excluding stars surrounded by extended Ha nebular 
emission. This prevents the removal of CVs with nova shells 
(for example T Aur in our In-Plane sample). 

4.2 Locating the Main Stellar Locus 

Objects surviving the initial cuts are then subjected to an it- 
erative selection algorithm (illustrated in Fig.0. First, they 
are divided into four magnitude bins, based on their r' band 
magnitudes. Next r' — Ha versus r' ~ i' colour-colour 
plots are created for each magnitude bin in each field. 

We next carry out an initial straight line least squares 
fit to the objects in each magnitude bin. Objects with strong 
Ha emission lines should exhibit an excess in r' — Ha and 
will therefore generally be located above the main stellar 
loci in the plots. The least squares fit is an initial attempt 
to fit the main stellar loci, and so we expect Ha emitters 
to reside above the fit in the plot. For fields with low num- 
ber density, which includes the majority of fields outside the 
galactic plane, this initial fit to all objects often works well. 
However, in fields with higher number densities, the stel- 
lar locus often splits due to the presence of different stellar 
populati ons (giants and dw arfs) and different reddening val- 
ues (see iDrew et al1l2005l for a detailed explanation). The 
key point for our purposes is that the upper locus of points 
in the colour-colour plots often runs above the initial least 
squares fit in fields dominated by more distant/highly red- 
dened population. Thus a selection criterion based solely on 
an object's position relative to our initial fit would result in 
the erroneous selection of many ordinary late-type unred- 
dened main-sequence stars which reside in the upper locus 
of points. 

4.3 Identifying The Upper Boundary of the Main 
Stellar Locus 

In order to prevent the selection of unreddened objects with 
no Ha emission, we use an iterative er-clipping technique to 
force the fit onto the upper boundary of the main stellar 
locus. The fact that the unreddened objects reside at the 
top of the main stellar loci is used in our technique to ex- 
clude gradually points below the initial fit line in subsequent 
iterations. In practice, we use four iterations to arrive at a 
final best-bet fit to the the upper boundary of the main stel- 
lar locus. The final iteration is used to give precedence to 
objects that are well separated from the initial fit (within 
reasonable limits) rather than those objects that lie in the 
regions where the upper boundary of the main stellar locus 
is no longer clearly separated from the lower regions of the 
main stellar locus. 

However, this final fit may not always be appropriate, 
such as in fields where the stellar locus is not actually split. 
In these cases, the gradient of the upper boundary of the 
main stellar locus fit often matches closely with the initial 
least squares fit or exhibits a nonphysically large slope. If 
either of these conditions is met, the selection of Ha emitters 
is based on the initial fit; in all other cases, the final fit is 
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Figure 1. An illustration of the selection criteria used to identify strong emission line objects on colour-colour diagrams. The data 
shown here are all from the IPHAS field 5455 which surrounds the bright dwarf nova EY Cyg. First, the data are split up into four 
magnitude bins, as shown in the four panels. Objects with Ha excess should be located near the top of a colour-colour diagram. The 
thick lines illustrate the fits to the data, and the thin lines illustrate the corresponding cuts used to single out clear Ha emitters. The 
solid lines show the original least squares fit to the full dataset, along with the corresponding selection cut. The dashed lines show the 
final fit and cut lines, obtained by applying an iterative cr-clipping techniques to the initial fit. The goal of the cr-clipping is to fit only the 
the upper boundary of the main stellar locus in the colour-colour diagram. Objects selected as Ha emitters are shown as large triangles. 
Note that the cut lines shown here are only approximate, as the actual selection criterion used also considers the errors on each individual 
datapoint (see Equation 121. 



used to select emitters. Empirically, we find that this method 
gives reasonable results 94 per cent of the time (compared 
to visual inspection) . That is, the fit runs through the upper 
boundary of the stellar locus when it is visible and otherwise 
fits the whole stellar locus in a satisfactory manner. 

4.4 Selecting Ha Excess Objects 

Once the appropriate fit for each bin has been decided, we 
identify objects significantly above the chosen fit as likely 
Ha emitters. Our selection criterion takes into account both 
the scatter of points around the stellar loci and the errors 
on the colours of each individual datapoint. We first define 
the Ha excess for each object as 

AHa = (/ - Ha) obs - (/ - Ha) flt (1) 



where (r'-Htt) b s is the observed value of the colour and 
(r'— Ha)fit is the value obtained from the fit for the corre- 
sponding observed value of r' — i' . The selection criterion is 
then given in terms of Ha excess by 

AHa > C v /K + ^. Ho) +ml t , ( a r , i , ) ]. (2) 

Here rms is the root mean square value of the residuals 
around the fit, C is a constant (we use values of 4.5 for our 
initial fits and 5 for the fits with cr-clipping), mat is the 
gradient of the fit line, and 07 r '-Ha) an d <Tr r '-i') are the 
errors on the observed colours. 

Note that the cut lines shown here are only approxi- 
mate, as the actual selection criterion used also considers 
the errors on each individual datapoint (see Equation 
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In the example shown in Fig.0 the lines illustrate the 
initial and final selection cuts applied. However the cut lines 
shown are only approximate, as the actual selection criterion 
used also considers the errors on each individual datapoint 
(see Equation and this is not reflected by the cut lines 
shown in the figure. This explains why in some cases objects 
appearing above the illustrated cuts are not selected as emit- 
ters. The selected Ha emitters are indicated by the large tri- 
angles. The figure shows data from IPHAS field 5455, which 
contains the CV EY Cyg. Note that the data quality of this 
field falls below the IPHAS standards, that is this field will 
be re-observed for the final IPHAS survey at a latter date. 
EY Cyg is a bright dwarf nova with V-band magnitudes of 
11.4 and 15.5 in outburst and quiescence, respectively. The 
IPHAS photometry (r 1 — 14.1) suggests the system was ob- 
served ne ar quiescence, when it is known to display clear Ha 
emission IIMunari et al.lll997f) . The C V is marked by a large 
open circle and was correctly selected as a likely emission 
line object by our algorithm. We note that, in this case, the 
final fit succeeded in identifying EY Cyg as an emission line 
object, whereas a selection based on the initial fit to all the 
data would have failed to do so. 

It is worth emphasizing that we make no claim that the 
method presented here is in any strict sense optimal for se- 
lecting Ha emitters. The design goals for our algorithm are 
merely that it should be simple, fully automated and pro- 
duce reasonable results in the vast majority of cases (when 
compared to a fully interactive outlier selection). Conse- 
quently, the selection criteria we use are essentially empirical 
and tailored to provide adequate results in most situations. 
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Table 1. Summary of the existing data for the CVs in the two samples. Magnitude ranges are drawn directly from the 
Downes et al catalogue. Orbital period information is drawn directly from the Ritter and Kolb catalogue. CV classification 
data are adapted from the information contained within CVCat. 



GOVS Name 


Other name 


Magnitude Range 


Orbital 


CV 


CV 










Period a 


Flag b 


Type c 










(days) 






Off-Plane 


RZ LMi 


PC 0948+344 


14.2 V 


- 16.8 V 


0.058500* 





DN/SU 


RU LMi 


CBS-119/Ton 1143 


13.8 p 


- 19.5 p 


0.251000 





DN 


CH UMa 


PC 1003+678 


10.7 V 


- 15.3 V 


0.343184 





DN 


SW Sex 


PG 1012-029 


14.8 B 


- 16.7 B 


0.134938 





NL/SW 


GG Leo 


RX J1015. 5+0904 


16.5 V 


- 18.8 V 


0.055471 





AM 


CP Dra 




14.3 p 


- 20 p 


0.081600* 





DN 


CI UMa 


SVS 1755 


13.8 p 


- 18.8 V 


0.060000 





DN/SU 


KS UMa 


RX J1020+5304 


13.0 p 


- 17.0 p 


0.068000 





DN 


FIRST J102347+003841 




17.3 V 


- 17.8 V 


0.198115 


-1 


- 


U Leo 


BD +14 2239 


10.5 v 


- <15 v 


0.267400: 





N: 


DW UMa 


PG 1030+590 


14.9 V 


- 18. V 


0.136607 





NL/VY 


DO Leo 


PG 1038+155 


16.0 B 


- 17.0 B 


0.234515 





NL 


IY UMa 


Tmz V85 


13.0 p 


- 18.4 p 


0.073909 





DN/SU 


CW 1045+525 




16.5 B 


- 


0.271278 





DN 


SX LMi 


CBS-31/Ton 45 


13 V - 


17.4 V 


0.067200 





DN/SU 


CY UMa 


SVS 2198 


12.3 V 


- 17.8 V 


0.067950 





DN/SU 


AN UMa 


PG 1101+453 


13.8 B 


- 20.2 B 


0.079753 





AM 


ST LMi 


CW 1103+254 


15.0 V 


- 17.2 V 


0.079089 





AM 


AR UMa 


1ES 1113+432 


13.3 V 


- 16.5 V 


0.080501 





AM 


DP Leo 


IE 1114+182 


17.5 B 


- 19.5 B 


0.062363 





AM 


TT Crt 


FSV 1132-11 


12.7 V 


- 16.3 V 


0.268420 





DN 


RZ Leo 




11.5 V 


- 19.2 V 


0.076038 





DN/SU 


T Leo 


BD +4 2506a 


10.0 V 


- 15.9 V 


0.058820 





DN/SU 


DO Dra 


YY Dra? 


10.6 B 


- 16.7 B 


0.165374 





IP 


TW Vir 


PG 1142-041 


12.1 v 


- 16.3 v 


0.182670 





DN 


EU UMa 


RE 1149+28 


16.5 B 


- 16.8 B 


0.062600 





AM 


BC UMa 


GR 102 


10.9 B 


- 18.3 B 


0.062610 





DN/SU 


IR Com 


S 10932 


13.5 B 


- 18. B 


0.087039 





DN:/SU: 


EV UMa 


RX J1307+53 


17.1 V 


- 20.8 V 


0.055338 





AM 


HV Vir 


NSV 6201 


11.5 V 


- 19.0 V 


0.057069 





DN/SU 


HS Vir 


PG 1341-079 


13.0 B 


- 16.6 V 


0.076900 





DN 


OU Vir 


1432-0033 


14.5 p 


- 18.5 j 


0.072730 





DN/SU 


UZ Boo 


HV 10426 


11.5 v 


- 20.4 V 


0.125000: 





DN/SU 


KV Dra 


RX J1450.5+6403 


11.8 V 


- 17.1 V 


0.058760 





DN/SU 


TT Boo 


HV 3681 


12.7 v 


- 19.2 V 


0.077000: 





DN/SU 


NY Ser 


NSV 6990 


14.8 V 


- 17.9 V 


0.097800 





DN/SU 


M5 V101 




17.5 p 


- 20.9 V 


0.242000 





DN 


ES Dra 


PG 1524+622 


13.9 p 


- 16.3 p 


0.179000 





DN/SU 


QW Ser 


Tmz V46 


12.8 p 


- <15.3 p 


0.074570 





DN 


ASAS 153616-0839.1 




11.5 V 


- <15.4 V 


0.063600* 





DN 


LX Ser 


Stepanian's star 


13.3 B 


- 17.4 B 


0.158432 





NL/SW 


CT Ser 


N Ser 1948 (r) 


7.9 v - 


16.6 p 


0.195000 





N 


SS UMi 


PG 1551+719 


12.6 V 


- 17.6 V 


0.067780 





DN/SU 


MR Ser 


PG 1550+191 


14.9 V 


- 17. V 


0.078798 





AM 



a A colon following the orbital period signifies the value is uncertain, an asterisk signifies the value has been estimated 
from the known superhump period. 

^ Possible entries in this column; = good, 1= uncertain, -1 = not in CV Cat, -2 = non cv. 

c Possible entries in this column; DN = dwarf nova, N = nova, NL = novalike variable, NR = recurrent nova, AM = 
polar (AM Her), IP = intermediate polar, - = non classified. Any entries in this column following a forward slash are the 
secondary classification of the CV, possible additional entries are; SU = SU UMa system, VY = VY Scl system, ZC = Z 
Cam system, UG = U Gem system, SW = SW Sex system, Na = fast nova, Nb = slow nova. A colon following any of these 
entries signifies the classification is uncertain. 
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Table 1 - continued 



GCVS Name 


Other name 


Magnitude Range 


Orbital 


CV 


CV 








Period a 


Flag b 


Type c 








(days) 






Off-Plane 


RX J1554.2+2721 




16.8 B - 


0.105462 





AM 


QZ Ser 


HadV4 


12.7 p - <14.9 p 


0.083161 





DN 


VW OB 


Antipin V21 


14.0 p - <17.5 b 


0.070700* 





DN 


1RXS J161008+035222 






0.132200 





AM 


X Scr 


HV 3137 


8.9 p - 18.3 p 


1.480000 





N:/Nb: 


V589 Her 


S 10296 


14.1 p - <17.5 p 


0.090500* 





DN 


V844 Her 


Antipin V43 


12.5 p - 17.5 p 


0.054643 





DN/SU 


V699 Oph 


HV 10577 


13.8 p - 18.5 p 


0.068300* 


-2 




RW UMi 


SVS 1359 


6 p - 18.5 B 


0.059100 





N/Nb 


V841 Oph 


BD -12 4633 


4.2 v - 13.5 v 


0.601400 





N/Nb 



Plane 



HT Cas 


S 3343 


10.8 v 


- 18.4 v 





073647 





DN/SU 


V Per 


BD +56 406a 


9.2 p - 


18.5 p 





107126 





N 


UV Per 


AN 87.1911 


11.7 V 


- 17.9 V 





.064900 





DN/SU 


TZ Per 


AN 28.1912 


12.3 v 


- 15.6 v 





262906 





DN/ZC 


T Aur 


BD +30 923a 


4.2 p - 


15.2 p 





.204378 





N/Nb 


FS Aur 


S 3946 


14.4 v 


- 16.2 v 





.059500 





DN 


CW Mon 


AN 61.1936 


11.9 v 


- 16.3 v 





.176600 





DN 


V603 Aql 


HD 174107 


-1.1 v - 


11.8 B 





.138500 





N/Na 


CI Aql 


AN 23.1925 


8.8 V - 


15.6 p 





.618363 





NR 


V446 Her 




3.0 p - 


17.8 B 





207000 





N/Na 


V1425 Aql 




7.5 v - 


<19.0 V 





.255800 





N/IP: 


EM Cyg 


AN 185.1928 


12.5 v 


- 14.4 v 





290909 





DN/ZC 


RX J1951. 7+3716 




15.5 B 







493000 







EY Cyg 


AN 200.1928 


11.4 v 


- 15.5 v 





.459326 





DN/UG 


V2306 Cyg 


1WGA J1958.2+3232 


15.7 V 




0. 


181200 





IP 


V550 Cyg 


S 3847 


14.2 p 


- <17.0 p 





.067200* 





DN/UG 


V503 Cyg 


S 4524 


13.4 V 


- 17.6 V 


0. 


.077700 





DN/SU 


V751 Cyg 


SVS 1202 


13.2 V 


- 16.0 V 


0. 


.144464 





NL/VY 


V2275 Cyg 


N Cyg 2001, No. 2 


6.6 p - 







.314500 





N 



a A colon following the orbital period signifies the value is uncertain, an asterisk signifies the value has been estimated 
from the known supcrhump period. 

^ Possible entries in this column; = good, 1= uncertain, -1 = not in CV Cat, -2 = non cv. 

c Possible entries in this column; DN = dwarf nova, N = nova, NL = novalike variable, NR = recurrent nova, AM = 
polar (AM Her), IP = intermediate polar, - = non classified. Any entries in this column following a forward slash are the 
secondary classification of the CV, possible additional entries are; SU = SU UMa system, VY = VY Scl system, ZC = Z 
Cam system, UG = U Gem system, SW = SW Sex system, Na = fast nova, Nb = slow nova. A colon following any of these 
entries signifies the classification is uncertain. 



Table 2: Summary of the IPHAS data for the CVs in the two samples. 
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Flag 


Flag 


(arcsec) 
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Off-Plane 


RZ LMi 


PG 0948+344 


09 


:51:48 


.93 


+34:07:23.9 


191.29 


51 


.04 


14 


789 


0. 


001 


14. 


783 





003 


14 


857 





.003 











1 


83 





09 


20 


.03 


-5 


RU LMi 


CBS-119/Ton 1143 


10 


02:07 


.47 


+33:51:00.5 


191.83 


53 


.17 


17 


626 


0. 


006 


17 


296 





013 


17 


007 





.008 


1 








1 


.42 





.10 


20 


.47 


65 


CH UMa 


PG 1003+678 


10 


07:00 


.59 


+67:32:47.2 


142.90 


42 


.65 


14 


327 





001 


13. 


683 





002 


13 


832 





002 


1 





1 


2 


.29 





.04 


20 


00 


35 


SW Sex 


PG 1012-029 


10 


:15:09 


.39 


-03:08:32.9 


245.52 


41 


.68 


15 


.071 


0. 


002 


14. 


952 





004 


14 


815 





.003 


1 





1 


2 


.53 





10 


19 


.62 


25 


GG Leo 


RX J1015. 5+0904 


10: 


15:34 


67 


+09:04:42.5 


231.59 


49 


.05 


18 


949 


0. 


014 


18 


936 





047 


19 


122 





.032 








1 


1 


.65 





.07 


20 


.20 


-15 


CP Dra 




10: 


15:39 


82 


+73:26:05.2 


136.34 


39 


.38 


15 


.587 


0. 


002 


15. 


628 





.005 


15 


796 





.004 











1 


.44 





07 


20 


.58 


-15 


CI UMa 


SVS 1755 


10: 


18:12 


.96 


+ 71:55:43.6 


137.59 


40 


.55 


18 


490 


0. 


010 


18 


247 





023 


18 


036 





014 


1 








1 


.20 





.12 


20 


.80 


45 


KS UMa 


RX J1020+5304 


10: 


20:26 


.52 


+53:04:33.2 


159.55 


51 


.93 


16 


.775 


0. 


004 


16. 


539 





008 


15 


953 





.004 
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.75 





09 


20 


.86 


110 


FIRST J102347+003841 




10: 


23:47 


.70 


+00:38:41.2 
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45 
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17 


315 
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16 


856 
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17 


315 
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1 


1 
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.07 


20 
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-10 


U Leo 


BD +14 2239 


10 


24:03 


92 


+ 14:00:25.3 


226.34 


53 


.27 


17 


007 


0. 


004 


16 


696 





009 


17 


063 





008 











1 


.08 





.10 


20 


.53 


-10 


DW UMa 


PG 1030+590 


10 


33:52 


86 


+58:46:54.8 


150.28 


50 


.41 


15 


211 


0. 


002 


15. 


116 





004 


14 


421 





.002 


1 











.70 





.07 


20 


.85 


110 


DO Leo 


PG 1038+155 


10 


40:51 


.23 


+ 15:11:34.0 


227.71 


57 


.43 


17 


591 


0. 


006 


16. 


947 





Old 


17 


136 





.009 


1 








1 


.32 





.04 


20 


.64 


30 


IY UMa 


Tmz V85 


10 


43:56 


72 


+58:07:31.9 


149.76 


51 


.84 


17 


638 


0. 


006 


17. 


299 





012 


17 


046 





.008 





2 


1 


2 


.31 





.12 


20 


.11 


65 


CW 1045+525 




10 


48:18 


01 


+52:18:30.0 


156.89 


55 


.92 


14 


608 


0. 


001 


14. 


000 





.002 


14 


164 





.002 


1 








1 


.50 





04 


20 


.63 


30 


SX LMi 


CBS-31/Ton 45 


10: 


: 54:30 


44 


+30:06:10.4 


199.34 


64 


.24 


16 


146 


0. 


003 


15. 


917 





006 


15 


596 





.004 


1 








1 


.36 





.05 


20 


.42 


60 


CY UMa 


SVS 2198 


10: 


56:57 


03 


+49:41:18.3 


159.37 


58 


.55 


17 


336 


0. 


005 


16. 


995 





011 


16 


534 





.006 


1 








1 


.39 





.05 


20 


.42 


100 


AN UMa 


PG 1101+453 


11 


04:25 


62 


+45:03:14.5 


165.83 


62 


15 


16 


.251 


0. 


003 


16. 


103 





006 


15 


908 





.004 


1 





1 


1 


.97 





.05 


20 


.33 


35 


ST LMi 


CW 1103+254 


11 


05:39 


79 


+25:06:28.7 


211.80 


66 


.21 


17 


523 


0. 


005 


16. 


681 





009 


17 


343 





.010 











1 


.30 





.07 


20 


.49 


-5 


AR. UMa 


1ES 1113+432 


11: 


:15:44 


56 


+42:58:22.6 


167.45 


64 


.97 


16 


.468 


0. 


003 


15 


962 





006 


16 


269 





.005 














.75 





.06 


20 


.77 


5 


DP Leo 


IE 1114+182 


11 


:17:15 


.93 


+ 17:57:41.8 


230.90 


66 


.46 


18 


082 


0. 


007 


17. 


915 





()18 


17 


723 





012 


1 





1 


2 


.06 





02 


20 


11 


35 


TT Crt 


FSV 1132-11 


11: 


34:47 


17 


-11:45:30.2 


274.88 


46 


.90 


15 


729 


0. 


002 


15. 


156 





004 


15 


379 





.003 


1 





1 


2 


.78 





11 


19 


48 


20 


RZ Leo 




11 


37:22 


.24 


+01:48:58.9 


264.77 


59 


.09 


18 


.715 


0. 


011 


18 


013 





.020 


17 


932 





.014 


1 





1 


2 


.47 





.03 


19 


.77 


80 


T Leo 


BD +4 2506a 


11: 


38:26 


81 


+03:22:07.5 


263.49 


60 


.52 


16 


048 


0. 


002 


15. 


676 





005 


14 


993 





.003 


1 





1 


2 


.69 





20 


19 


.42 


155 


DO Dra 


YY Dra? 


11 


43:38 


46 


+ 71:41:20.9 


130.30 


44 


.46 


14 


516 


0. 


001 


13. 


872 





002 


13 


485 





.001 


1 








1 


.37 





07 


20 


.53 


130 


TW Vir 


PG 1142-041 


11: 


45:21 


.18 


-04:26:05.7 


273.60 


54 


.63 


15 


763 


0. 


002 


15. 


199 





004 


15 


112 





.003 


1 








1 


.74 





.05 


20 


.33 


60 


EU UMa 


RE 1149+28 


11: 


49:55 


.71 


+28:45:07.5 


202.54 


76 


33 


17 


283 


0. 


004 


17. 


140 





011 


16 


477 





.006 


1 





1 


2 


.49 





05 


19 


.87 


110 


BC UMa 


GR 102 


11 


52:15 


.84 


+49:14:41.9 


146.27 


65 


.12 


18 


157 


0. 


007 


17. 


971 





019 


17 


570 





.011 


1 





1 


2 


.04 





11 


20 


.01 


70 


IR Com 


S 10932 


12 


39:32 


.05 


+21:08:06.3 


277.93 


83 


.42 


18 


618 


0. 


010 


17. 


809 





018 


18 


534 





.020 











1 


.08 





.04 


20 


55 


-10 


EV UMa 


RX J1307+53 


13 


07:53 


82 


+53:51:31.0 


117.57 


63 


.10 


16 


630 





003 


16. 


759 





009 


16 


495 





.006 


1 








1 


.07 





06 


20 


.65 


20 


HV Vir 


NSV 6201 


13: 


21:03 


.18 


+01:53:29.0 


319.88 


63 


.78 


19 


058 


0. 


014 


19 


013 





040 


18 


938 





025 





2 





1 


.18 





.07 


20 


74 


10 


HS Vir 


PG 1341-079 


13 


43:38 


44 


-08:14:03.6 


324.44 


52 


.44 


15 


206 


0. 


002 


14. 


980 





003 


14 


839 





.003 


1 








1 


.17 





12 


20 


.74 


35 


OU Vir 


1432-0033 


14 


35:00 


.23 


-00:46:06.2 


348.90 


52 


.60 


16 


556 


0. 


003 


16. 


445 





007 


16 


013 





.005 


1 





1 


1 


.84 





.06 


20 


.11 


65 


KV Dra 


RX J1450. 5+6403 


14 


50:38 


34 


+64:03:28.8 


103.81 


48 


.41 


17 


086 


0. 


004 


16. 


766 





009 


16 


439 





.006 


1 








1 


.80 





03 


20 


.25 


70 


TT Boo 


HV 3681 


14 


:57:44 


79 


+40:43:40.9 


68.74 


60 


.70 


13 


688 


0. 


001 


13. 


722 





002 


13 


740 





.002 





2 








.74 





.07 


20 


.99 


-5 


NY Ser 


NSV 6990 


15 


13:02 


.33 


+23:15:08.5 


33.99 


57 


.84 


15 


311 


0. 


002 


15. 


269 





004 


15 


341 





.003 














.82 





.04 


20 


.82 


-5 


ES Dra 


PG 1524+622 


15 


:25:31 


.80 


+62:00:59.8 


97.64 


46 


83 


15 


.019 


0. 


001 


14. 


823 





003 


14 


905 





.003 


1 











.89 





17 


20 


.72 


5 


QW Ser 


Tmz V46 


15 


26:13 


.99 


+08:18:02.5 


13.06 


48 


.86 


17 


603 


0. 


005 


17. 


140 





.012 


16 


941 





.008 


1 








1 


.72 





.04 


20 


08 


70 


ASAS 153616-0839.1 




15 


36:16 


.00 


-08:39:07.7 


356.94 


36 


.40 


18 


.096 


0. 


007 


17 


852 





018 


16 


993 





.008 


1 











.80 





.11 


20 


.76 


175 


LX Ser 


Stcpanian's star 


15 


38:00 


11 
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29.78 


50 
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14 
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14 
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-5 
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15 
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.35 


+27:21:52.7 


44.16 


49 


.61 


17 
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16 
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16 
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17 
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0. 


005 


16. 


262 





007 


16 
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31 
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17 
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16 
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30 
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16 
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17 


610 


0. 
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2 
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16 


:25:01 
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17 
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16 
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Table 2 - Continued 
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-2.71 


15 


936 





004 


15 


455 





005 


15 


056 
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01 
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17 
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133.47 
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4.5 Interpreting Ha Excess and Estimating 
Equivalent Widths 

It is tempting to treat the Ha excess values calculated from 
Equation 1 as direct tracers of emission line strength (trac- 
ers of line EW). However, this interpretation is not always 
safe. In particular, it is important to keep in mind that the 
calculated excess values have no fixed reference point, that is 
each excess is calculated relative to the local fit to the upper 
boundary of the stellar locus in the appropriate magnitude 
bin of a given field. Thus any direct comparison of Ha ex- 
cess values between different objects implicitly assumes that 
each object actually belonged to the stellar population that 
was described by the fit. In fields containing a split stellar 
locus (and hence more than one population), there is no 
guarantee that this is the case. 

It is worth stressing that this problem is not simply due 
to the particular definition of "Ha excess" we have adopted. 
Instead, it reflects the inherent ambiguity in assigning indi- 
vidual objects to particular stellar populations on the basis 
of limited photometric information. More specifically, since 
different stellar populations can have different "baseline" 
(r'-Ha) colours at given (r' — i'), it is impossible to de- 
fine a physically meaningful Ha excess for a given source 
without specifying the population to which it belongs. 

This all sounds rather pessimistic. However, we can also 
turn this argument around: if we are willing to adopt a par- 
ticular spectral energy distribution (SED) for a given Ha 
excess source, a corresponding Ha EW estimate can be ob- 
tained immediately from the photometric data. If this SED 
represents the objects'continuum distribution better than 
the fit to the upper stellar locus, the resulting EW estimate 
should be a more meaningful quantity than the straight Ha 
excess. 

Fig.|5| illustrates how this idea can be implemented for 
CVs. Basically, the location of a given source in the colour- 
colour diagram is determined by its Ha EW and the shape 
of its continuum SED. The latter, in turn, depends on the 
intrinsic SED shape and the reddening towards the source. 
For CVs, it may be reasonable to adopt an intrinsic SED of 
Fx cc A" 2,3 as typical. This is appropriate for an optically 
thick accretion disc and is probably acceptable even for WD- 
dominated short-period DNe. We can then add Ha lines of 
varying EWs to this SED and calculate the reddening tracks 
of the resulting model spectra through colour-colour space. 
As shown in Fig. [5] this procedure yields a mapping between 
the observed colours, on the one hand, and the correspond- 
ing Eb-v values and Ha EWs, on the other. 

In principle, it should be possible to place observed 
points directly onto this mapped colour-colour plot and read 
off the resulting Eb-v values and Ha EWs. However, in 
practice, one additional correction is necessary at present. 
This arises because our synthetic photometry is Vega-based, 
but the standard stars used for the photometric calibration 
of the IPHAS data a re mostly later type stars. As noted 
bv iDrew et alJ (|2005l) . Vega is an AOV star with deep Ha 
absorption, and this causes a systematic offset between ob- 
served and synthetic r' - Ha colours. The offset required for 
a given field is between 0.25 and -0,05 mag (in the sense 
that the data points must be shifted upwards). We thus ac- 
count for this by applying an average shift of 0.15 mag to 



Table 3. The distribution of Ha emitters amongst the CV sam- 
ples. The numbers are given as a ratio of CVs selected as Ha 
emitters to the total number of CVs in the sample. The numbers 
in parenthesis show the corresponding percentages. 



Sample 


All CVs 


CVs Passing 






initial cut 


Off-Plane 


36/52 (69.2) 


36/48 (75.0) 


In-Plane 


12/19 (63.2) 


12/17 (70.6) 


Total 


48/71 (67.6) 


48/65 (73.8) 



the observed r' - Ha colours of all CVs before placing them 
on Fig.H 

It is also worth noting that the IPHAS data cur- 
rently still lack a global (field-to-field) photometric calibra- 
tion. This again particularly affects the Ha photometry and 
causes a scatter of ~ 0.1 mag in the r' - Ha colours between 
different fields. We currently make no attempt to correct for 
this. Note that both systematic and field-to-field offsets will 
be automatically corrected in the final IPHAS data release. 

Fig-El shows the corrected colours of our CVs on the 
mapped colour-colour diagram. The resulting Ha EW esti- 
mates are given in Table El and we will adopt these values 
throughout the rest of this paper. We stress that these may 
be subject to considerable error, both because of field-to- 
field scatter and because the intrinsic SED shape we have 
assumed may not be appropriate to all systems. 4 



5 RESULTS 

Table summarizes the key properties of all the CVs in the 
two samples. The magnitude range data were tak en from 
the Downes et al catalogue 5 jDownes et al.ll200ll) . orbital 
periods are from the Ritter & Kolb catalogue, and the CV 
type informati on comes fr om the web-based CV catalogue 
CVCat 6 jKube et al.l2003l) . Our new data are listed Tabled 
which also includes information on whether a given CV was 
selected as an Ha emitter. 

5.1 Recovery Rates of CVs as Ha Emitters 

Our recovery rates are summarized in Tables 13151 and illus- 
trated in Figs 13141 In the following subsections, we discuss 
the resulting completeness estimates and their dependence 
on key CV properties (orbital period, CV type, apparent 
magnitude, galactic latitude). We also explore the causes 
that prevent some CVs from being selected as Ha emitters. 

5.1.1 Overall Completeness 

Table El gives the recovery rates (that is CVs selected as 
Ha emitters) for both in-plane and off-plane samples and 
also the combined recovery rates. We find a global recovery 

4 On a more positive note, Figure 6 in lDrew et alJ J2005I) shows 
that, for smooth SEDs, most of the sensitivity to the assumed 
SED shape is in the reddening estimates, whereas Ha EWs are 
much less affected. 

5 http:/ /archive. stsci.edu/prepds/cvcat/index. html 

6 http://cvcat.net 
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Figure 2. Colour-colour plot showing how Ha EW relates to Ha excess and broad- band colour. Triangles represent dwarf novae; circles 
represent novae, novalike variables and recurrent novae; squares represent magnetic systems and crosses represent non-classified systems. 
The symbols are plotted in red if the system is a Ha emitter; blue if it is a non-emitter; and green if it is a system which did not pass 
our initial selection cuts. The jagged solid curve is the unreddened main sequence (MS) line, the smooth solid curve running mainly 
horizontally is the early-A reddening line (that is the reddening track of a n early A star). A ll other lines correspond to power law SEDs 
with F\ oc A~ 2 3 , as appropriate for an optically thick accretion disc (see lDrew et alJl2005l) . The dashed lines are lines of constant Ha 
EW superposed on this SED (the EW values are given in Angstroms at the far left of each line). The solid green vertical lines represent 
Ha emission curves of growth for this SED at fixed reddening (the corresponding E(B - V) values are shown at the top of each line). 



rate of ~ 68% [48/71] for all CVs and ~ 74% [48/65] if we 
only consider those CVs that passed our initial selection cuts 
defined in section 14.11 

Both of these numbers are important. The first pro- 
vides an "end-to-end" estimate of the survey completeness 
for CVs, that is what fraction of CVs brighter than r' — 19.5 
should we expect to discover with IPHAS (or similar sur- 
veys) as emission line objects. This estimate allows for in- 
completeness due to crowding/blending, objects falling on 
bad pixels, astrometric calibration problems, etc. By con- 
trast, the second completeness estimate (considering only 
objects that pass our initial selection cuts) is an estimate 
of the fraction of CVs that are in principle discoverable by 
IPHAS, that is the fraction of CVs that exhibit photometri- 
cally detectable Ha excesses. Note that the number of CVs 
detected as Ha emitters is the same in both cases, since a 
CV which does not pass the initial selection cuts can never 
be selected as an Ha emitter. 

It is interesting to ask what constitutes a "detectable" 
Ha excess in this context. A useful perspective on this ques- 
tion can be obtained from Fig.|5| which shows that the 



unreddened main sequence broadly separates sources de- 
tected as Ha emitters from sources not detected as Ha excess 
objects. This reflects the fact that our selection procedure 
uses the upper boundary of the main stellar locus as the zero 
point in searching for "significant" Ha excess outliers and 
confirms that this boundary usually tracks the unreddened 
main sequence. 



It is worth noting that this effectively constitutes an ob- 
servational selection bias. For a given SED, reddened sources 
must have larger Ha EWs than unreddened sources in or- 
der to be selected as Ha excess objects. Most CVs are blue 
objects observed through relatively little extinction, so the 
sense of this selection bias works in their favour. In essence, 
CVs occupy a region in the colour-colour plots where even a 
small Ha EW stands out prominently from the stellar locus. 
The high recovery rate we find is partly due to this fortunate 
situation. 
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Figure 3. a) The magnitude distribution of the CVs b) the period distribution of the CVs. The orbital period gap is shown between 
the cyan lines. The plots show the distribution of emitters in red overlaid on top of the total distribution of CVs. 



Table 4. A summary of the results showing the distribution 
amongst the long and short period CVs. Numbers are given as a 
ratio of CVs selected as Ha emitters to the total number of CVs 
in the sample, the percentages are given parenthesis. 





Sample 


All CVs 


CVs passing 








initial cut 


Short 


Off-Plane 


21/34 (61.8) 


21/31 (67.7) 


Period 


In-Plane 


4/5 (80.0) 


4/5 (80.0) 


< 2.5 hr 


Total 


25/39 (64.1) 


25/36 (69.4) 


Long 


Off-Plane 


15/18 (83.3) 


15/17 (88.2) 


Period 


In-Plane 


8/14 (57.1) 


8/12 (66.7) 


> 2.5 hr 


Total 


23/32 (71.9) 


23/29 (79.3) 



5.1.2 Dependence on Orbital Period 

In order to assess the dependence of our completeness on 
orbital period, we first compare the recovery statistics of 
short-period CVs (defined as having orbital periods below 
the 2.2 - 2,8 hour period gap) to those of long-period CVs 
(with periods above the gap) . These recovery rates for both 
types of CVs are listed in Table [I] 

Somewhat surprisingly, the recovery rates for long- and 
short-period systems are comparable, with both lying near 
65 per cent. Given that short-period CVs tend to be intrin- 
sically faint and that the H/3 EW and absolute magnitude 
for DNe are correlated (in the sense that the faint systems 
have stronger lines; Patterson 1984), one might have ex- 
pected that short-period systems would be easier to detect 
and thus have higher recovery rates. The fact that we do 
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Table 5. Results showing the distribution of Ho emitters 
amongst the two samples of CVs. The results in this case are 
separated by CV type. Numbers are given as a ratio of CVs se- 
lected as Ho emitters to the total number of CVs in the sample, 
and the percentages are shown in parenthesis. 



CV type Sample All CVs CVs passing 

initial cut 

Dwarf Off-Plane 22/30 (73.3) 22/27 (81.5) 

Novae In-Plane 6/9 (66.7) 6/9 (66.7) 

Total 28/39 (71.8) 28/36 (77.8) 

Magnetic Off-Plane 7/11 (63.6) 7/11 (63.6) 

Systems In-Plane 1/1 (100.0) 1/1 (100.0) 

Total 8/12 (66.7) 8/12 (66.7) 

Nova, Recurrent Off-Plane 7/9 (77.8) 7/8 (87.5) 

Nova and Novalike In-Plane 4/8 (50.0) 4/6 (66.7) 

Variables Total 11/17 (64.7) 11/14 (78.6) 

Non Off-Plane 0/2 (0.0) 0/2 (0.0) 

Classified In-Plane 1/1 (100.0) 1/1 (100.0) 

Total 1/3 (33.3) 1/3 (33.3) 



not observe a marked difference in recovery rates is mainly 
a testament to the high data quality of IPHAS, that is its 
ability to separate even long-period objects with relatively 
weak emission lines from the main stellar loci in any given 
field. This explanation is supported by the fact that short- 
period DNe do exhibit the largest Ha EWs (see section 16.11 
and particularly Fig.|SJ. 

5.1.3 Dependence on CV Type 

In order to examine whether Ha selection favours particu- 
lar CV sub-types, each system in our samples is assigned 
one of four basic types based on their classification in CV- 
Cat shown in Tabled (i) dwarf novae; (ii) novae and nova- 
likes; (iii) magnetic CVs; (iv) unclassified systems. Note that 
the second class includes both classical and recurrent no- 
vae, and the third class includes both intermediate polars 
and the strongly magnetic polars. We refrain from using a 
finer classification scheme since this would result in classes 
with extremely small numbers. It should also be pointed out 
that some intermediate polars are known to undergo dwarf 
nova outbursts. For the purposes of the work presented here, 
any system which has been classified as a magnetic CV is 
put into this category, irrespective of any similarities with 
weakly-magnetic systems. 

Table |^| shows the recovery statistics as a function of 
CV types and Fig. [3] shows how the recovery statistics de- 
pend on CV magnitude and orbital period for the different 
CV types. The dwarf novae have the overall highest recov- 
ery rate, which is not surprising since, in quiescence, the ac- 
cretion discs in these systems are at least partially optically 
thin and are w ell known to produce strong emission lines (see 
Warner 1995 and references therein). The magnetic systems 
show the second highest recovery rate of all systems. This is 
probably due to strong Balmer emission from their accretio n 
streams dCropperlll99tilF'errario fc Wickramasinghelll993l) . 

The novae and nova-likes have the lowest overall recov- 
ery rates, and it is tempting to attribute this to their known 
observational and physical characteristics. At first sight, this 



Table 6. A summary of the results, in this case showing how 
the recovery rate of Ho emitters varies between the four magni- 
tude bins. The numbers are given as a ratio of CVs selected as 
Ho emitters to the total number of CVs in the sample and the 
percentages are shown in parenthesis 



Mag 
rang 


nitude 

;e (r') 


S<im.plc 


All CVs 


CVs passing 
initial cut 


< 16 


Off-Plane 

In-Plane 

Total 


10/16 (62.5) 
7/11 (63.6) 
17/27 (63.0) 


10/14 (71.4) 
7/10 (70.0) 
17/24 (70.8) 


16 - 


17.5 


Off-Plane 

In-Plane 

Total 


14/19 (73.7) 
2/3 (66.7) 
16/22 (72.7) 


14/19 (73.7) 
2/3 (66.7) 
16/22 (72.7) 


17.5 


- 18.5 


Off-Plane 

In-Plane 

Total 


10/12 (83.3) 
3/5 (60.0) 
13/17 (76.5) 


10/11 (90.9) 
3/4 (75.0) 
13/15 (86.7) 


18.5 


- 19.5 


Off-Plane 

In-Plane 

Total 


2/5 (40.0) 
0/0 (N.A.) 
2/5 (40.0) 


2/4 (50.0) 
0/0 (N.A.) 
2/4 (50.0) 



seems reasonable: both types of high-state, weakly-magnetic 
CVs tend to exhibit only weak emission lines when viewed 
at moderate inclinations (for example IShara et al.l Il98fil : 
IWarnerlll987l : lDhillonlll99d) . This is probably because these 
systems contain optically thick discs that produce mainly 
absorption line spectra (although often emission reversals 
are observed in the cores of absorption line profiles). 

However, a closer look at the statistics suggests that the 
true explanation for the low recovery rates for these systems 
is more prosaic. In the last column of Table we show the 
recovery rates that result if we restrict attention to systems 
that have passed the initial selection cut. Here, the recov- 
ery fraction of novae and nova- likes is the highest. These 
numbers are completely consistent with the corresponding 
ones for DNe and magnetic systems. Thus the reason for the 
lower overall recovery statistics for novae and nova-likes is 
simply that relatively more of these systems are misclassified 
as non-stellar in the IPHAS data. 

It is worth asking why, despite the physical arguments 
given above, the recovery rates for novae and nova-likes are 
comparable to those of DNe. This question has actually al- 
ready been answered implicitly in section 15.1.21 There, we 
noted that short-period CVs (which are dominated by DNe) 
do exhibit larger Ha excesses than long-period CVs (which 
contain most of the novae and nova- likes), but that IPHAS 
is sensitive enough to still identify the weaker Hcv emission 
often seen in the latter. 

Finally, the non-classified systems systems show a com- 
paratively low recovery rate in comparison to other systems. 
Given that there are only 3 non- classified CVs it is not pru- 
dent to draw any firm conclusions from the low recovery 
rate. 

5.1.4 Dependence on CV Magnitude 

As photometric scatter increases towards fainter magni- 
tudes, completeness may be expected to decrease. Clearly, a 
faint CV will on average need to display a greater Ha excess 
in order to be detected. Table|S|shows how the recovery rate 



14 A.R. Witham et al. 



of Ha emitters depends on the magnitudes bins used in the 
selection. It can be seen that for the Off-Plane sample the 
recovery rate actually increases slightly as we go to fainter 
bins until we get to the faintest bin, where it drops to 40 
per cent (but note that there are only 5 sources in this bin). 
The In-Plane sample also hovers at a near constant 60 to 
70 per cent for the three brighter bins, and for this sample 
the faintest bin contains no objects. We conclude that com- 
pleteness for CVs in IPHAS remains steady and high up to 
at least r' ~ 18.5. The small numbers of fainter CVs in our 
samples prevent us from constraining the drop in complete- 
ness beyond this point more precisely. 

5.2 Reasons for CVs Not Being Selected as Ha 
Emitters 

There are several reasons that can prevent a CV from being 
selected as an Ha emitter: 

(i) The CV is not located significantly above the main 
stellar locus in the colour-colour plots. In some cases, this 
may be due to weakness or absence of an Ha emission line (as 
seen in erupting dwarf novae and some nova-like variables). 
Furthermore, the equivalent width of the Ha emission line 
can sometimes change significantly during the orbital cycle, 
so some CVs may have been observed at phases where the 
emission is weak. Finally, at fainter magnitudes, the signal- 
to-noise may simply become insufficient to separate a CV 
from the main stellar locus. 

(ii) The CV is located below the locus in the colour-colour 
plots, that is the CV shows Ha absorption. This may again 
point towards an erupting dwarf nova or a nova-like variable. 

(iii) The fitting process failed, and the CV is not selected 
due to the resulting selection cut being inappropriate. 

(iv) The CV does not pass the initial selection cuts de- 
fined in section 14.11 

In the following sections, we look in more detail at the 
various CV sub-types and check if we can understand why 
particular CVs were not recovered as emission line objects. 

5. 2. 1 Dwarf Novae 

The DNe are the most numerous class and generally show 
some of the largest Ha excesses. However, they also show the 
greatest range in Ha excesses, although this can be at least 
partly attributed to observing some systems in quiescence 
and others in outburst. RZ LMi, CP Dra, TT Boo, NY Ser, 
and TZ Per are the DNe that have been observed in out- 
burst. None of them have been detected as Ha emitters. Of 
the 6 DNe remaining that were not selected as Ha emitters, 
IY UMa has a large Ha excess but does not pass the initial 
selection cuts, for the reasons given in section 13.11 and hence 
is not selected. CW Mon also has a large Ha excess but is 
not selected as described in section 16.11 HV Vir misses our 
initial selection cuts and hence is not selected, however its 
EW shows that Ha emission is present. This leaves 3 systems 
which are not recovered as Ha emitters: IR Com, VW CrB 
and V550 Cyg. These systems pass our initial selection cuts 
and do not have a field quality flag associated with them. 
It is likely that variations in EW due to orbital phase can 
account for the absence of observable Ha emission in these 



systems, for example IR Com is known to be an eclipsing 
system and the IPHAS magnitudes are consistent with this 
system being observed in eclipse. 

The recovery rate for all the DNe in quiescence is 81 
per cent 7 which is 9 per cent higher than the recovery rate 
for all the DNe including those in outburst. 

5.2.2 Magnetic CVs 

There are rather few magnetic systems in our samples (in- 
cluding only one intermediate polar), but the majority show 
a detectable Ha excess. Those which do not show emission 
are all AM Her systems observed in the low state (GG Leo, 
ST LMi, AR UMa, and MR Ser). A reduction in Ha emission 
line strength is a known characteristic of low state magnetic 
CVs. It is possible that the Ha emission is not strong enough 
to be detected above Zeeman absorption in these systems. 

5.2.3 Novae, Recurrent Novae and Novalike Variables 

Many of the bluest CVs of this class show a large Ha ex- 
cess indicating that the nova remnants are undergoing mass 
transfer. However, there is no simple link between detection 
of Ha excess and time since eruption. For example, V Per 
and T Aur erupted in 1891 and 1887, respectively, and both 
have been detected as emitters. Conversely, the two oldest 
nova remnants - V841 Oph (1848) and U Leo (1855, though 
its nature as a nova is not certain) - have not been detected 
as Ha excess objects. 

There is one nova with a very high Ha excess that did 
not pass the initial selection cut. This system is V1425 Aql. 
Spectroscopic observations hav e shown strong Ha e mission 
both before and after outburst llKamath et al.lll997f) . which 
is consistent with the large excess we observe. 

5.3 In-Plane vs Off-Plane CVs 

Given that IPHAS is a galactic plane survey, but that the 
majority of the CVs analyzed here are Off-Plane sources, it 
is worth checking whether crowding and extinction affects 
the recovery rate in galactic plane fields. We have already 
provided separate recovery rates for both samples in Ta- 
bles but in Fig.^Jwe now additionally present a visual 
comparison between the recovery statistics of In-Plane and 
Off-Plane samples. 

In general, we find that the recovery rates of the two 
samples are similar. This suggests that crowding and ex- 
tinction have at most a small effect on completeness. How- 
ever, despite the similar recovery rates, there is nevertheless 
a clear difference in data between galactic plane fields and 
non-galactic plane fields. For example, when the value of the 
rms scatter about the initial least square fits to the stellar 
loci is compared between In-Plane and Off-Plane fields, the 
scatter is clearly larger for the former. This is mainly due 
to different populations seen at different distances through 
varying extinction in galactic plane fields (though probably 

7 The CVs HS Vir and OU Vir have been excluded from this 
recovery rate calculation as the photometry suggests they are in 
an intermediate state between quiescence and outburst. 
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Figure 4. The magnitude distribution and period distribution of CVs showing those located within the galactic plane (upper plots) 
and outside of the plane (lower plots). The plots show the distribution of emitters in red overlaid on the total distribution of CVs. The 
orbital period gap is shown between the cyan lines. 



also to increased photometric scatter due to higher crowd- 
ing levels) . Our cr-clipping fits to the upper boundary of the 
main stellar locus help partly to counter this effect. With- 
out this technique, the In-Plane recovery rates would be 
reduced. 



6 CORRELATIONS BETWEEN PHYSICAL 
AND PHOTOMETRIC PROPERTIES 

6.1 Orbital Period vs Ha EW 

In Fig. |^| we plot the Ha EW - estimated as described in 
Section 14.51 - versus the orbital period for each CV. As ex- 
plained in the introduction, short-period CVs are expected 
to show stronger Ha emission, and this is indeed confirmed 
by the figure. 

The colour codes and plot symbols are defined in the 
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Figure 5. A plot of EW versus the orbital period of CVs in the Off-Plane and In-Plane samples. The orbital period gap is shown 
between the cyan lines. 



figure legend and are used to distinguish between the CV 
types (symbols) and the emitters (colours). A small dot in 
the centre of each symbol refers to the data in the CV flag 
and CV Type column of table The dot is coloured ma- 
genta if the CV flag is not equal to zero (the CV is not in 
the catalogue or it is of uncertain nature) or the CV type 
is of an uncertain nature. If the CV flag is zero and the CV 
type is certain, the dot is black. 

The short-period quiescent DNe show the greatest EW, 
and the survey has no problems in detecting them as emis- 
sion line objects. In fact the long-period nova V1425 Aql 



shows the biggest EW but it misses our initial selection cuts 
and the EW estimate is too large to be physical, hence it is 
not shown in the figure. There are more short-period DNe 
than long-period DNe in our sample, which is a characteris- 
tic of the Ritter & Kolb catalogue as a whole (and expected 
on evolutionary grounds). 

As already noted in section 15.1.21 long-period systems 
have a lower Ha excess, but the survey is still sensitive 
enough to detect many of them as emitters. As shown in 
Fig-El the boundary between systems detected as emitters 
and those undetected corresponds to a Ha EW of ~ 10A; 
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the average excess for all systems is clearly above this value, 
irrespective of period. Thus most long-period CVs do ex- 
hibit photometrically detectable Ha emission, albeit at a 
lower EW level than short-period systems. 

There is one object in Fig.[I)]that shows a large apparent 
Ha EW, but was not selected as an emitter. This object is 
the long-period dwarf nova CW Mon, which is located in a 
field of poor data quality. In fact, this field does not meet 
the IPHAS criteria for inclusion in the final survey and will 
be re-observed. As a result, the stellar locus exhibits a large 
scatter which places the cut further up in the colour-colour 
plots compared to fields with good data quality. The CV is 
separated from the locus to the extent that, on our colour- 
colour plots, it appears above the graphical representations 
of the selection cuts. However, it is not selected due to the 
large uncertainties of the magnitudes. Past spectroscopic ob- 
servations hav e revealed strong Balmer emission in CW Mon 
l|Szkodvlll987h . 

6.2 r' - i' vs Orbital Period 

In Fig.^we plot r' — i! colour against orbital period. No ob- 
vious correlation is seen, but we note that the r' — i' colours 
of the CVs in this plot span a range of more than 1 magni- 
tude. Based on this, we believe the use of selection criteria 
based on this broad-band colour (such as a "blue cut" ) would 
not be prudent in the construction of a new CV sample from 
IPHAS. In this context, it is also worth keeping in mind that 
the sample of known CVs analyzed here may already suffer 
from selection effects. It is therefore important to minimize 
the number of a priori selection criteria when constructing 
new samples that are meant to better represent the true CV 
parent population. 

6.3 Ha EW vs Absolute Magnitude 

In order to test for a correlation between intrinsic bright- 
ness and EW, we calculated absolute magnitudes for CVs 
with known distances in our sample by assuming a galactic 
extinction law of Ay = 1.5 mag/kpc. In doing so, we also 
assumed an Ry =3.1 extinction curve to convert Ay to A r i . 
Fig.[7| shows the relationship between the resulting absolute 
magnitudes and our estimated EWs for various groupings. 

In the plot showing all systems, there appears to be a 
correlation, albeit with large scatter. The existence of a cor- 
relation becomes more apparent when we plot the individual 
CV types separately. For the DNe and the novae, in partic- 
ular, there is a definite trend for the Ha EW to increase 
towards larger (that is fainter) absolute magnitudes. The 
magnetic systems do not show this trend, which is reason- 
able given the lack of accretion discs in the AM Her systems. 

Our r esults agree qualitatively with the work of 
iPattersonl (ll984T> . who has calculated a fit to the relation- 
ship between absolute magnitude and H/3 equivalent widths. 
Note, however, that Patterson removed the contribution of 
the secondary from his absolute magnitude estimates and 
thus considered only the intrinsic brightness of the accretion 
disc in his correlation. We have not made this correction, nor 
have we corrected in any way for inclination. Nevertheless, 
all of our results clearly confirm that short-period, faint CVs 
(which are expected to represent the most numerous CV 



population) exhibit the strongest Ha emission lines. Thus 
photometric emission line surveys are indeed an excellent 
resource for finding and studying this population. 



7 DISCUSSION 

7.1 Prospects for CV Searches Based on Emission 
Line Surveys 

A key question we wanted to address with this study was 
"can surveys like IPHAS find members of the postulated 
large population of faint, short-period CVs"? The results 
here confirm that short-period CVs generally exhibit rela- 
tively large Ha excesses in IPHAS, so if the hypothesized 
large population of short-period CVs conform to this rule, 
the survey would have no problem detecting them. More sur- 
prisingly, we have also been able to recover long-period CVs 
as Ha emitters with a comparable success rate. Thus we ac- 
tually do not expect any additional bias from the survey in 
favour of detecting either long- or short-period CVs. 

The vast majority of short-period systems are believed 
to be DNe. However, due to the low mass accretion rates in 
these systems, the outburst duty cycle is extremely low. For 
example, the proto-type WZ Sge erupts on average every 
30 yrs and stays bright for only a few months. In the con- 
text of an emission line survey for CVs, this may actually be 
an advantage. The results here show that IPHAS is unable 
to detect most DNe in outburst via our selection method, 
whereas quiescent DNe tend to show the greatest Ha ex- 
cesses. Thus short-period DNe are actually more likely to 
be discovered by IPHAS if they have long outburst inter- 
vals. 

It is interesting to ask how many faint, short-period CVs 
IPHAS might eventually discover. Given the low accretion 
rates and very low mass secondaries, it is possible that the 
optical luminosity of these systems will actually be domi- 
nated by their WD primaries. An accretion-heated WD in 
a short-period systems would be expe cted to have an effec- 
tive t emperature of T e ff ~ 10000 K llTownslev fc Bildstenl 
2003). IPHAS could see such an object out to 275 pc, 
given its limiting magnitude of r' = 19.5. For an assumed 
CV space density of 5 x 10 -5 pc -3 (see discussion below), 
there are roughly ~ 0.1 short-period CVs per square degree 
brighter than r' = 19.5 mags near the galactic plane. Since 
IPHAS will image 1800 square degrees, we would predict 
that ~ 180 of these systems could be detected with IPHAS. 

This prediction is conservative in the sense that no ac- 
tive CV should be fainter than the estimate used above - 
any admixture of brighter CVs will increase the total num- 
ber considerably (there are intrinsically fewer of these, but 
they can be seen to greater distances). The prediction does, 
of course, scale directly with the assumed space density. 
The number adopted above is based o n theoretical popu- 
lation synthesis predictions dKolbl 19931). and is higher than 
most existing empirical estimates Patterson 1998; G ansickd 
2005). This is not a problem, as it is precisely this discrep- 
ancy between theory and observation that an IPHAS CV 
sample would be designed to test. 

It has already been noted in Section I5.1.4l that a faint 
CV will on average need to display a greater Ha excess than 
a bright CV in order to be detected. However, if the trend 
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Figure 6. A plot of r' — i' versus the orbital period of CVs in the Off-Plane and In-Plane samples. The orbital period gap is shown 
between the cyan lines. 



towards larger EWs for short-period, faint DNe continues to 
hold for this population, then there may be no problem in 
detecting these systems even close to the magnitude limit of 
the survey. 



Considering all of the above, it seems clear that IPHAS 
(and similar surveys) does have the potential to test theo- 
retical population synthesis predictions. 



7.2 Implementing an IPHAS-based CV Search 

Initial follow-up to the IPHAS survey has already lead to 
many clear Ha emitters being spectroscopically observed 
and identified. A key problem in constructing a new sample 
of CVs in this way is the presence of other large populations 
of objects displaying Ha emission. Selecting objects for spec- 
troscopic follow-up using the selection techniques illustrated 
here can uncover most of the CVs, but the fraction of CVs 
amongst objects selected in this way may not be large. Our 
results to date suggest that early type emission line stars 
(for example Be stars) tend to dominate bright (r' ^ 18) 
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Figure 7. The relationship between the absolute magnitude and the Ha equivalent width for the CVs with known distances. The 
symbols are plotted in red if the systems is a Ha emitter; blue if it is a non-emitter; and green if it is a system which did not pass our 
initial selection cuts. 



samples of Ha excess objects. Thus many non-CVs need to 
be observed in the process of constructing a new sample of 
CVs, unless the distribution of Ha excess objects changes 
significantly towards fainter magnitudes. We are currently 
investigating this possibility. 

8 CONCLUSIONS 

We have considered the properties of known CVs contained 
within the galactic plane Ha survey IPHAS as of July 2004, 



along with those of a sample of CVs outside the galac- 
tic plane observed with the same observational setup. We 
have developed routines to select clear Ha emitters from 
the IPHAS photometry and used these routines to deter- 
mine the percentage of CVs in our samples that would be 
recovered as Ha emitters. 

We find that the overall completeness of IPHAS-like 
surveys for detecting CVs as emission line objects is ~ 70 
per cent, roughly independent of CV type and orbital pe- 
riod. Our recovery fractions are steady near this value up to 
at least r' ~ 18.5 and only drop towards the very faintest 
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magnitudes we consider (r' = 19.5). There are 23 CVs in 
our samples not recovered as Ha emitters. Six do not pass 
our initial selection cuts, and 2 are without a known sub- 
type. Several are not recovered due to observing some CVs 
when Ha emission is expected to be weak or absent, for ex- 
ample DNe in outburst, and magnetic CVs in a low state. 
The remainder of non-recovered CVs are due to scatter in 
faint magnitude bins and cases where the fitting process did 
not produce good results. 

We have estimated Ha EWs from our photometry and 
find that, as expected, short-period CVs show the strongest 
Ha emission lines. This is encouraging for the prospects of 
finding the predicted population of faint, short-period CVs. 
However, we also find that most long-period CVs exhibit 
strong enough Ha emission for them to be easily detected 
by IPHAS-like surveys. 

In closing, we expect that IPHAS (and related surveys) 
will allow us to construct samples of CVs that are free from 
most of the selection effects that have plagued previous com- 
parisons between theoretical population synthesis predic- 
tions and observations. In particular, such samples should 
uncover the long- predicted large population of faint, short- 
period systems, if it exists. 
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